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Abstract Angular momentum in protostellar discs can 
be transported either radially, through turbulence in- 
duced by the magnetorotational instability (MRI), or 
vertically, through the torque exerted by a large-scale 
magnetic field. We present a model of steady-state discs 
where these two mechanisms operate at the same ra- 
dius and derive approximate criteria for their occur- 
rence in an ambipolar diffusion dominated disc. We 
obtain "weak field" solutions - which we associate with 
the MRI channel modes in a stratified disc - and trans- 
form them into accretion solutions with predominantly 
radial angular-momentum transport by implementing a 
turbulent-stress prescription based on published results 
of numerical simulations. We also analyze "intermedi- 
ate field strength" solutions in which both radial and 
vertical transport operate at the same radial location. 
Our results suggest, however, that this overlap is un- 
likely to occur in real discs. 

Keywords accretion, accretion discs — MHD — stars: 
formation — ISM: jets and outflows. 

1 Introduction 

Angular momentum in protostellar discs can be trans- 
ported either radially, through turbulent viscosity in- 
duced by the magnetorotational instability (MRI; e.g. 
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Balbus & Hawley 1998); or vertically, by centrifugally 
driven winds (e.g. Konigl & Pudritz 2000). It is likely, 
in fact, that both forms of transport play a role in 
real discs Quasi-steady disc models that incorporate 
the two mechanisms have already been attempted (e.g. 
Lovelace et al. 1994; Casse & Ferreira 2000; Ogilvie 
& Livio 2001), but these studies treated the radial 
transport using prescriptions that do not account for 
its origin in MRI-induced turbulence. A more refined 
approach seems to be needed, however, as solutions 
in which a wind carries away all the excess angular 
momentum of the fluid are not susceptible to sustain 
MRI-unstable modes (Wardle & Konigl 1993, hereafter 
WK93). Our prescription to incorporate radial trans- 
port in our wind-driving models uses the fact that the 
growth (and structure) of the MRI is dependent on the 
local strength of the magnetic field, which we measure 
by a, the ratio of the Alfven speed to the isothermal 
sound speed (see section |2|. MRI modes are expected 
to grow when a ^ 1 (provided that the gas is suffi- 
ciently coupled to the field), but they are suppressed 
when a is comparatively large (> 1). Since a (oc p""'^) 
increases with height, this suggest that radial and ver- 
tical angular momentum transport could in principle 
operate at the same radius, with the MRI acting in a 
finite section next to the midplane, and outflows dom- 
inating at higher z. 

Since protostellar discs are weakly ionized over most 
of their extent, it is essential to incorporate the ef- 
fect of finite magnetic diffusivity in our models. This 
is done in our general scheme by employing a tensor- 
conductivity, obtained via a realistic vertical ionization 
structure. This formulation naturally incorporates the 



^ Note that vertical angular-momentum transport could also oc- 
cur via other mechanisms, not considered here; such us magnetic 
braking (e.g. Krasnopolsky & Konigl 2002), "failed" winds or 
non-steady phenomena. 
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three relevant magnetic diffusion meclianisms (i.e. am- 
bipolar, flail and Ohniic). The illustrative examples 
presented here, however, assume that ambipolar diffu- 
sion dominates over the entire section of the disc and 
(except in section |4]) , take the magnetic coupling to be 
constant with height. This is useful in order to apply 
the analytic results obtained by WK93 for this case to 
our prescription for MRf-induced turbulent transport 
in a wind-driving disc model. We refer the reader to 
Salmeron, Konigl & Wardle 2007 (hereafter SKW07) 
for full details of the study. Section [2] presents the 
wind-driving model. Section [3] outlines our prescrip- 
tion for MRI-induced radial transport and illustrates 
it with a solution with a predominantly radial angular- 
momentum transport. Section|4]discusses a representa- 
tive solution that incorporates radial and vertical trans- 
port and our findings are summarised in section [sj 

2 Centrifugally driven winds 

Bipolar outflows are commonly associated with proto- 
stellar accretion discs, as well as other astronomical ac- 
creting sources. Their ubiquity suggests that they play 
a key role in regulating accretion, a notion supported by 
the apparent correlation between accretion and outflow 
signatures in these systems (e.g. Hartigan, Edwards & 
Ghandour 1995). It is also generally accepted that these 
outflows are highly coUimated winds, accelerated cen- 
trifugally from the surfaces of discs by the mechanism 
first proposed by Blandford & Payne (1982, hereafter 
BP82). According to this picture, disc material is flung 
out if the open magnetic field lines that thread the disc 
are sufficiently inclined (i.e. the angle between the line 
and the disc plane is < 60°), initiating an outffow that 
can potentially reach super Alfvenic speeds. 

A key feature of these winds, and one of particu- 
lar interest in the context of this work, is that they 
provide an efficient means of extracting the rotational 
kinetic energy of the disc and carrying away the ex- 
cess angular momentum of the accreted matter. This is 
expected to be the dominant angular momentum trans- 
port mechanism in discs when the magnetic field is 
strong (suprathermal) and the MRI is suppressed be- 
cause the wavelength of the most unstable perturba- 
tions exceed the magnetically reduced disc scale height. 
Wind-driving solutions from strongly magnetised discs 
have been obtained semi-analytically and numerically 
(e.g. WK93; Li 1995, 1996; Casse & Keppens 2002). 

In this work we have used the modelling procedure 
developed by WK93 to obtain radially locaHsed disc- 
wind solutions from strongly magnetised discs. This 
entails solving the conservation (of mass and momen- 
tum) and induction equations, as well as Ampere's and 



Ohm's laws for a steady-state, isothermal, geometri- 
cally thin and nearly Keplerian accretion disc. We ne- 
glect all radial derivatives except that of and use a 
conductivity tensor (e.g. Wardle 1999) to communicate 
the effect of the charged species on the neutrals. 

These equations are integrated vertically upwards to 
obtain the disc solution and the location of the criti- 
cal (sonic) point. This solution is then matched onto 
a BP82, self-similar wind solution, by imposing the 
regularity condition at the Alfven critical point. The 
relevant model parameters are (1) oq = wao/cs, the 
midplane (subscript 0) ratio of the Alfven speed to 
the sound speed, which measures the magnetic field 
strength; (2) rj, the ratio of the Keplerian rotation time 
to the neutral-ion momentum-exchange time, which de- 
termines the magnetic coupling; (3) e = —Vro/cs, the 
midplane (normalized) inward radial speed, evaluated 
by imposing the Alfven critical-point constraint on the 
wind; (4) Cs/vk = ^T/r, the ratio of the disc tidal scale 
height to the radius; and (5) es = —cE^/csBz, the 
normalized azimuthal component of the electric field, 
which is nonzero if the magnetic field lines drift radi- 
ally (WK93)|3 

The main properties of wind-driving, strongly mag- 
netized (ao ^ 1) discs are analysed in WK93. Fig. 1 of 
SKW07 shows a solution of this type that matches on 
to a global BP82 wind solution (see also the pure-wind 
case of Fig. [2]). These solutions satisfy: 

{27j)-^/^ < ao < V3 < er] ; rj > 1 . (1) 

The first inequality in ([T]) expresses the requirement 
that < vk in the disc region, the second gives the 
wind launching condition {B^ i,! > l/\/3, where sub- 
script 'b' denotes the disc surface), the third constrains 
the base of the wind to lie above the magnetically- 
reduced density scale height {h) and the fourth is the 
minimum field-matter coupling condition. 

3 MRI-induced turbulence 

The MRI transports angular momentum radially out- 
wards as fluid elements exchange angular momentum 
non-locally via the distortion of the magnetic field lines 
that connect them. Its properties have been exten- 
sively investigated in both the linear and nonlinear 
stages (e.g. Sano & Stone 2002a,b; Salmeron & War- 
dle 2003, 2005). In particular, numerical simulations 



^Setting £3=0 effectively fixes tlie value of at the disc sur- 
face (subscript 'b'). In a global formulation Br would be deter- 
mined by the radial distribution of B^ (e.g. Ogilvie & Livio 2001; 
Krasnopolsky & Konigl 2002). 



Radial and vertical angular momentum transport in protostellar discs 



3 



show that the MRI sustains radial angular-momentum 
transport when the initial (subscript i) Elsasser num- 
ber Aj = v\q -Jriohm^K is 1 (e.g. Sano & Inutsuka 
2001; Sano & Stone 2002a, b), where "qohm is the Ohmic 
diffusivity. For the vertical initial field geometry as- 
sumed here, it turns out that A; = rji (SKW07) and 
the minimum-coupling condition is the same for both 
forms of angular- momentum transport. So the ques- 
tion we seek to answer is: Assuming that i] > I, which 
section of the disc (if any) would be unstable to the 
MRI? The first inequality in ([ij can in fact be used 



m 
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Fig. 1. — Structure of a weakly magnetized disc with 
vertical magnetic angular-momentum transport in the 
limit where no wind develops. The radial and vertical 
components of the magnetic field are indicated by the 
thick solid and dashed lines, respectively. The vertical 
dotted line marks the boundary of the MRI-unstable 
zone (see text). 

to derive a lower limit of a for MRI stability. This is 
the case because when this expression is violated, the 
surface layers of the disc become super-Keplerian and 
experience outward streaming motion that is unphysi- 
cal for a pure wind solution (WK93), but is character- 
istic of the two-channel mode that underlies the MRI. 
Moreover, wind-driving disc solutions that satisfy all 
the constraints in ([iJ are stable to the fastest grow- 
ing linear mode of the MRI (WK93), a condition that 
further supports this choice. The above constraint on a 
was applied by WK93 at the midplane, since in their so- 
lutions the wind carries away the entire excess angular 
momentum of the gas. Here, however, it is more appro- 
priate to apply it locally [using a{z) cx p(z)~^/^] to dif- 
ferentiate the MRI-unstable section of the disc (where 
2ria^ < 1) from the region where only vertical angular 
momentum transport takes place {2rjo? > 1). 

Fig. [l] shows a solution where a < (2ry)~^/^ over 
most of the disc vertical extent. Note that no wind 



develops and thus, angular momentum is transported 
magnetically between different heights and ultimately 
carried away radially by the fiow. As anticipated, both 
the radial and azimuthal field components oscillate in 
the region where 2r]a^ < 1 and the disc sections where 
BzdB^/ dz < (> 0) lose (gain) angular momentum 
and have < (> 0). We associate this solution 
with a two-channel MRI mode for a stratified, non-ideal 
MHD disc. These modes were analysed by Goodman & 
Xu (1994) in the unstratified, ideal-MHD case. These 
authors showed that they are exact solutions in the lin- 
ear and nonlinear regimes, but are unstable to parasitic 
modes. Although a rigorous stability analysis of the 
modes in stratified, non-ideal MHD discs have not yet 
been done, numerical simulations suggest that they may 
indeed evolve into turbulence (e.g. Hawley et al. 1995; 
Sano et al. 2004, hereafter S04). 

Having derived a criterion to identify the MRI- 
unstable section of the disc {rj > \ and 2ria? < 1), we 
now need to quantify the turbulent angular momentum 
transport that we expect to develop in these regions. To 
accomplish this, we use the results of numerical sim- 
ulations in the literature to evaluate the spatial and 
temporal average of the r<j) component of the turbulent 
Maxwell stress (the main component of the stress ten- 
sor) in terms of the similarly averaged magnetic pres- 
sure, obtaining 



< Wr0 > ~ 0.5 



< 5^ > 
Stt 



(2) 



(e.g. Hawley et al. 1995; S04). Using the characteris- 
tic ratios for the field components in MRI turbulence 
reported in S04, this can be expressed as 



< Wrd, >~ 14 F 



Bl 



Stt 



(3) 



where F ee< > /B^ ^ jg also listed in S04 as a 
function of /3i = 2/a?, the initial plasma /3 parameter. 
Finally, the radial angular-momentum transport associ- 
ated with this term is added to the angular momentum 
conservation equation in the disc region deemed to be 
unstable to the MRI. This equation then reads 



PVrVK 

2r 



dz c r 



(4) 



where we used (l/r)9(r^ <C Wr^, ^)/dr ft!<C 
d{rv^) jdr k, vk/2 and « 0. Note that we are implic- 
itly assuming that the turbulent magnetic field compo- 
nents only contribute to the radial angular-momentum 
transport and interpret the other magnetic-field (and 
velocity) terms in the disc structure equations as the 
mean values of the respective quantities. Similarly, we 
take the first term on the r.h.s. of equation (|4| as the 
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vertical angular-momcntuni transport induced by the 
mean field. The underlying assumptions of this ap- 
proach would need to be verified numerically. 

4 Combined angular momentum transport 

We now explore solutions where both radial and vertical 
angular momentum transport take place at the same 
radial location in the disc. From the above discussion, it 
is clear that they require moderate coupling (rj ~ 1) and 
field strengths [ao < (2ryo)~^/^ but a{z = zi,) ^ 1] for a 
wind to develop and the MRI to be active over a sizable 
section of the disc. To satisfy these requirements, we 
have allowed rj to decrease with z and fall slightly below 
1 in these solutions (with the expectation that it is still 
large enough to sustain MRI turbulence; see Sano & 
Stone 2002a, b). The adopted profile (with 77 decreasing 
from 0.63 to 0.5 between the midplane and the surface) 
is motivated by the inferred behavior of ri(z) in the 
ambipolar-diffusion dominated regions of protostellar 
disks (see SW07 and Salmeron & Wardle 2005). 

Fig. [2] shows solutions obtained under these assump- 
tions. The solid lines display a pure wind (verti- 
cal angular-momentum transport only) solution. The 
dashed lines depict how this solution is modified when 
radial angular-momentum transport is incorporated via 
the prescription introduced in Section 3. Both solutions 
are compared in Table [l] Note that the addition of ra- 
dial angular-momentum transport results in this case in 
an increased inflow speed (measured by e) and a higher 
mass inflow rate (Afjn). The higher \vr\ corresponds 
to a stronger radial neutral-ion drag (for cb = 0', see 
footnote |2| , and therefore to a stronger radial bending 
of the magnetic field lines (e.g. a higher B^.b)- This, in 
turn, increases the level of magnetic squeezing (reflected 
in the lower value of h), and leads to a stronger density 
stratification (e.g. a lower Ps/po)- Although the disc 
thickness does not change appreciably (zb = 1.25 /it in 
both cases) and the height of the sonic point (zs/Ht) 
decreases (from 2.90 to 2.82), the lower Ps/po leads to a 
lower wind outflow rate and vertical torque (measured 
by K and Tz, respectively) in the combined solution. 

The above results suggest that the joint operation 
of these two angular-momentum transport mechanisms 
might be self-limiting (in which case it would be un- 
likely to introduce a dynamical instability): The in- 
corporation of MRI-induced turbulence increases the 
overall angular momentum transport in the disc, but 
reduce the amount of angular momentum carried away 
by the wind. This conclusion could be tested through 
an explicit stability analysis (e.g. Konigl 2004) or via 
numerical simulations. 



Table 1 Comparison of the solutions in Fig. [2] Mjn 
is the mass accretion rate per disc circumference, nor- 
malized by (pqCs/it)- The vertical (Tz) and radial (Tr) 
torques per unit area are normalized by roBQ/An. All 
heights are listed in units of h = z{p = pi^j \fe). 



Disc / outflow 


Pure wind 


Combined 


properties 


solution 


solution 




0.31 


0.29 


Zh/h 


3.99 


4.35 


Zs/h 


9.25 


9.81 


PsIPa 


1.0 X 10-2 


8.3 X 10-3 


K 


0.17 


0.14 


Br.h/Bz 


1.75 


1.81 


\B(ji^b\/Bz 


1.44 


1.43 


Tz 


1.63 


1.59 


Tr 




0.09 


A/in 


1.74 


1.84 



5 Discussion 

We have examined the two leading processes thought 
to be responsible for transporting angular momentum 
in protostellar discs: Jets and outflows accelerated cen- 
trifugally from the disc surfaces and turbulent viscosity 
induced by the magnetorotational instability. Our main 
motivation has been twofold. First, we have set out 
to explore these processes under more realistic condi- 
tions, taking into account the detailed ionisation struc- 
ture and vertical stratification of the disc. Second, we 
are interested in the possibility that both forms of trans- 
port could operate at the same radial location and, if so, 
what is their relative importance at different locations 
and evolutionary stages in protostellar discs. Our over- 
all formulation incorporates a height dependent, tensor 
conductivity. This enables us to model the different 
diffusion mechanisms available to the fluid (ambipolar, 
hall and ohmic) and how their relative contributions 
change with location. The illustrative examples pre- 
sented here, however, assume ambipolar diffusion dom- 
inates over the entire section of the disc and the mag- 
netic coupling is constant with height. Our main results 
are highlighted below. 

1. We derived an approximate criterion to identify 
the section of the disc that would be unstable to 
MRI-induced turbulence and developed a scheme, 
based on the results of published numerical sim- 
ulations, to quantify the amount of angular mo- 
mentum transported radially in this region. 

2. We showed a weak field solution exhibiting oscil- 
latory field and velocity profiles and associated 
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Fig. 2. — Wind-driving disc with a moderately strong field and r] decreasing from 0.63 to 0.5 between z = and 
z = Zf,. Density, magnetic field components (left) and velocity components (right) are shown for a pure wind 
solution (solid lines) and also for a case where radial transport is included (dashed lines) in the region where 
2r]a'^ < 1, or to the left of the vertical dotted lines. The curves terminate at the respective sonic points. 



it with the nonlinear MRI channel mode identi- 
fied by Goodman & Xu (1994) in an ideal-MHD, 
unstratified case. A rigorous stability analysis of 
our sohitions has not been attempted yet, but ev- 
idence from linear analyses and numerical simu- 
lations suggests that they may be unstable. 

3. We argued that radial angular-momentum trans- 
port operates where 2?7a^ < 1 and presented a hy- 
brid solution in which this transition occurs at an 
intermediate height below the launching region. 

4. This solution shows that the addition of radial 
transport increases the total amount of angular 
momentum removed from the inflowing gas, but 
reduces the angular momentum flux carried away 
by the wind. 

5. It appears that these two transport mechanisms 
are unlikely to have a significant spatial overlap in 
real discs, as this requires a moderate «o {}t 0.5) 
and low (and decreasing with height) values of 
V ivo ^ 1) that are unlikely to be attained over 
measurable radial extents. 

We are currently in the process of applying our gen- 
eral scheme to the investigation of various conductiv- 
ity regimes in realistic disc models. This would en- 
able us to estimate the fractions of angular momen- 
tum transported via turbulent viscosity and outflows, 
respectively; and assess the implications of the result- 
ing disc structure to other discs processes, in particular 
to planet formation and migration. 
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